clc;
clear;
fh=@(x)sqrt(1-x.^2/2);
fl=@(x)-sqrt(1-x.^2/2);
f=@(x,y)exp(-x.^2/2).*sin(x.^2+y);
y=integral2(f,-1/2,1,fl,fh)


syms x y
s=int(exp(-x^2/2)*sin(x^2+y),y,-sqrt(1-x^2/2),sqrt(1-x^2/2));
int(s,x,-1/2,1),
vpa(ans)